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Abstract 

Examples of scalar conserved currents are presented for trigonometric, hyper- 
bolic and elliptic versions of the Hubbard model with non-nearest neighbour variable 
range hopping. They support for the first time the hypothesis about the integra- 
bility of the elliptic version. The two-electron wave functions are constructed in an 
explicit form. 
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The Hubbard models with non-nearest-neighbour hopping first proposed by Gebhard 
and Ruckenstein [1-4] have received considerable attention since this family of models 
itself contains the usual Hubbard model on an infinite lattice as well as the Haldane- 
Shastry [5,6] spin chain as certain limiting cases, and so do the generators of its Yangian 
symmetry [7]. Here we would like to present some conserved currents belonging to the 
scalar part of the algebra which provide an evidence for the integrability of more general 
models with hopping given by elliptic functions. In paper [7] it has been shown that 
trigonometric and hyperbolic models with hopping matrices tjk = t(j — k) = sin -1 jr(j — k) 
and tjk = t(j — k) — sinh^ 1 n(j — k) have an explicit Yangian symmetry which, however, 
does not imply integrability. One needs also the scalar part of the whole symmetry 
algebra which has been found in [8] for the Haldane-Shastry spin chain and in [9] for the 
conventional Hubbard model. 

The basic model describes electrons of spin a = ±1/2 created by operators Cj a at site 
j of a one dimensional lattice. The probability amplitude for hopping between sites j 
and k will be denoted by tjk- The strength of the repulsive interaction of two electrons of 
different spin on the same lattice site is U > 0. In these notations, the Hamiltonian reads 

h = £*i*<&Q* + 2C/£(4t c ;t - ±)(c+ Cji - J). (i) 
j,k j 

We shall consider translational invariant hopping amplitudes, tjk = tj-k- 
There is a canonical transformation, namely 

C jl -> C 3h C it C jP U ~ U - ( 2 ) 

This transformation leaves every Hamiltonian of the form (1) with antisymmetric hop- 
ping matrix invariant, but the global spin operators and the Yangian generators are not 
invariant under (2). 

To write down the SU2 generators of the rotational symmetry of the Hamiltonian ([I]) 
and generators of other symmetries, it is convenient to use spin operators formed as linear 
combinations of products of one creation and one annihilation operator at one site. It 
is also useful to introduce spin-like operators with indices corresponding to two different 
sites so as one arranges the pair of operators Cj a Ck T in a 2 x 2-matrix labeled by spin 
indices a and r, {Sjk)% = Cj a Ck T , 

S°<k = tr(a* a Sjk), S% = tr(Sjk), S? := S%, fij = Sj,, (3) 
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where the a a are the Pauli matrices. Note that (5" fc ) + = S k j, (Sj k ) + = S k j. The spin 
density and electron density operators are denoted by and S®, respectively. The 
commutators of these operators can be written explicitly, 

[SjkiSlm] = $kiSj m - S mj Si k , (4) 
[SjkiSfm] = 3klSf m — 5 m jS? k , (5) 

[Sf k ,sL] = 5^(5 kl S^ m -5 mj Sf k )+is a ^(5 kl S] m + 5 mj Sl). (6) 

However, there are other bilinear relations due to the composite nature of these operators. 
The Hamiltonian ([[]) now takes the following form 

^ = E^ + ^E((^-i) 2 -l)- (7) 

j,k j 

Since the particle number 7° = J2j Sj is conserved, only the term (Sj) 2 is relevant in the 
interaction part of the Hamiltonian while the other terms commute with H and can be 
removed by a shift of the chemical potential. 

To provide examples of the conserved currents which might exist for some choice of 
the hopping matrix, consider the ansatz 

N 



J = E [ A Jk S jk + B jk(SjS k - SjS k ) + D jk (Sj + S k )Sj k + E jk (S? v ' 



in which N is the number of sites. The condition [H, J] = with the use of (4-6) can be 
cast into the form of two functional equations 

4t jk (B lk - Bji) + (t jt D lk - Djitik) = 0, (9) 
2(t jk E kl + t kj E jl ) + (t jl D kl + t kl D jl ) = 0. (10) 

The coefficients A can be expressed in terms of B, D and E as 

Aj k = —2Dj k + (2U) [—8tj k Bj k + 2t k jEj k — aj k ), (11) 

where 
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Several boundary equations for t, B and D must be satisfied, too: 

= 0, 
= 0, 
tkjAjk) = 0. 



(tjiAik — Ajitik) 
^2(tjkD kj - D jk t k j) 



The first functional equation is just the Calogero-Moser one with known general analytic 
solution if;: 

a{x + A) 
nX) a(x)a(\f ■ 

By using it, one can express all the unknown structures in (9-11) (recall that tjk = t(j — k) 
etc.) as follows, 

t(x) = t ip(x), B(x) = --if;(x)ip(-x), 



D{x) = d 



, E{x) 



dip 2 (x) 



1 - hip(x+\)if>(-x-\) 



a(x) = t d^(x) — (N — 3)p(x) + hi(N - 2)£(x) + (£(x) - h 1 )(2xC(N/2) - N((x)) + s 



N-l 



e(x) = C(x + A)-C(x)-C(A), fci=V(A)/2, s = -(N-2) p (X)-J2p(1), 



i=i 



where a, ( and p are the Weierstrass elliptic functions determined by the full periods oj\ = 
N, uj 2 = in- The other parameters are given by A = ia or ia + N/2, v — — 2((N/2)\N~ 1 , 
in which k, d, h and a are arbitrary real numbers. Besides this general solution, there are 
of course the degenerate rational, trigonometric and hyperbolic forms which correspond 
to one or two infinite periods of the Weierstrass function. The vanishing of the boundary 
terms can be verified rather easily for the degenerate cases; in general elliptic case one 
needs to perform long but straightforward calculations which show that the boundary 
terms also disappear at the above choices. The key formula in these calculations is the 
following identity 

[p(y + A) - p(\)\ [C(x -y)- ((x + A) + ((y) + C(A)] + 

[p(x + A) - p(X)} [C(y -x)- C(y + A) + ((x) + C(A)] = p'{\). 

One thus can see that the formulas for t, B, D, E and a define scalar conserved currents 
for two three- parametric families of Hubbard models defined by the sets (A = ia, k, U/t ) 
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and (A = N/2 + ia,K,U/to). At the points A = wi/2,cj 2 /2 or (u\ + uj 2 )/2, i.e. at A 
being half-periods of the Weierstrass function p, the function ip(x) becomes odd and 
another independent current appears under the canonical transformation (2). It would be 
of interest to verify by direct calculations that both currents commute as it takes place 
for the two copies of the Yangian operators in the case of trigonometric and hyperbolic 
hopping [7]. One can thus conclude that these three-parametric families of Hubbard 
models might be integrable. An important open question is to confirm this conjecture 
by constructing the complete set of conserved currents commuting with the Hamiltonian. 
One can also see that the trigonometric and hyperbolic models of Bares, Gebhard and 
Ruckenstein fall into these families if one of the periods of the Weierstrass function, uj\ or 
u>2, tends to oo under an appropriate choice of A. 

Let us now calculate exact two-electron wave functions for the Hamiltonian (1) with 
the general elliptic hopping matrix. The problem is nontrivial only for the 5 = sector 
of the model. The wave functions can be written as 

* = £^+4ilo> + £<^+ T c+|o>, (12) 

where |0) is the vacuum state. The eigenequation HQ = can be written in the form 

£ (tjs<f>sk + tks<f>js) + {tjk9k + tkj9j) = E<f) jk , (13) 

£*;*(&* + = (E-2U) 9j . (14) 

Mi 

Note that 0^ can always be decomposed into a sum of the symmetric and antisymmetric 
parts 

fak = ^(<Pjk + <Pkj) + -((/)jk ~ <Pkj)- 

For antisymmetric <f>jk, the equation (14) would be satisfied trivially by gj = 0, and the 
solution to the equation (13) would be given by antisymmetrized product of one-electron 
plane waves. Thus from now on we concentrate on the nontrivial symmetric part, or in 
other words, assume that 4>jk is symmetric: 

4>jk = 4>kj- 

The ansatz for (f>j k and gj reads 

jfc = e'friJ+w*)^ -k) + e i( P lk+P2j ^ip(k - j), (15) 
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9j = 9oe 



i(pi+P2)j 



where 



<PU) = V?o 



(16) 



The parameters p 1 and p 2 are connected with r by the conditions of the periodicity of 

4>jk, 4>j+N,k = 4>j,N+k = 4>jk, 



(17) 



where rji = 2((N/2). The problem is to find all the parameters pi j2 , go and r from the 
eigenequations (13-14) if the ansatz (15-16) is correct. Recall that the elliptic hopping 
matrix in general is given by 



t(j — k)— t ij}{j — k)— t e 



uu-k) a U -k + X) 



(18) 



a(j-k)a(\y 

where the factor v = —2((N/2)\N~ 1 is chosen so as to satisfy the periodicity condition, 
t(j — k + N) = t(j — k). With the use of (16) and (18) the second eigenequation (14) can 
be cast into the form 



(E - 2U)g = 2to<p Si(pi,P2,T), 



(19) 



where 



S 1 (p 1 ,p 2 ,r) = XX 



a(s + X) 



-IP2S 



a(s + T) + e' tpiS a(s-T 



7#> ff2 ( s V( A ) 
The first eigenequation can be written as 

F(r,p 1 ,p 2 ,j,k) + F(r,p 1 ,p 2 ,k,j) + F(-r,p 2 ,Pi,j,k) + F(-r,p 2 ,Pi,k,j) (20) 



= -go 



t(j - k)e l{pi+P2)k + t(k - j)e 



i(pi+P2)j 



+ E<fi jk: 



where 



F(r, Pl ,p 2 ,j,k)= £ t(j-s)e^ s+p ^ V (s~k) 
= e i(pij+P2k) t(-q)e ipiq cp(q + j-k). 
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Let us now calculate the sum in the last expression with the use of the explicit forms of 
(16), (18), and introduce the notation 

N-l 



[21) 



where I = k — j G Z, and the function of a continuous argument x, 

G(l, x) = £ a{q - X \ X ]^ q 7 l+ /\ x) e^K (22) 

The function G(l,x) is double quasiperiodic, 

G(l,x + 1) =e u - tpi G{l,x), G{l,x + u 2 ) = e" 2(T_A) G(Z, x), 

where 771 = 2((N/2) and 772 = 2£(cj2/2) (recall that c^2 is the second period of all the 
Weierstrass functions here). It has a simple pole at x = with decomposition near it of 
the form 

Gfti)|„„= - (*-' + CH + r)- CH) - C(A)) 

On the other hand, this function can be written in the form 

G(i,x) = G (iy* ai(x + fl \ 

(Ti{X) 

where o\ is the Weierstrass sigma function with quasiperiods (1, 0U2) and Gq(1) is a constant 
factor. The parameters r and ji can be found from the quasiperoidicity conditions, 

r(r,pi) = (2m)~ 1 (r] 12 {u - ipi) - (r - A^^n) , 
V>(t,Pi) = (2m)- 1 (- W2(v - ipi) + (r- A)r/ 2 ) , 

where 7711 = 2^(1/2) and 77 12 = 2( 1 (u; 2 / / 2). Comparing the decompositions of both forms 
near x — 0, one finds the sum <S(Z) explicitly, 

m = - a{ ~!V ] (cm + r - a-i + r) + c(-o + c(a)) (23) 



(t(-J) 
<t(Z - A)ct(t) 

^5Ra) 



(Ci(/i) + r - c(i - A) + c(0 - COO) 



The explicit form of F(r,pi,p 2 , j, k) now reads with the use of equations (22), (23) as 
follows, 



F(r,p 1 ,p 2 ,j, k) 




Now, taking explicit summation of all four F's with different arguments in equation (20), 
one finds its compact form 



where r = r(— r, p 2 ), ft, = fi(—r,p 2 ). Comparing the coefficients in both sides, one obtains 
two equations 



These equations define E and <7o/Vo m terms of p±,p 2 and r. Plugging them into (19) 
results in the equation for determining the phase shift r. Together with the relations 
between r,p\ and p 2 (17) this equation allows one to determine all the parameters of the 
ansatz (14) for two-electron wave functions. 

To summarize, we have found new scalar conserved currents for the Hubbard model 
with an elliptic hopping matrix and its trigonometric and hyperbolic degenerations. We 
proposed the ansatz for the two-electron S = wave function and proved that it allows 
one to determine momenta and phase shift of two-electron states. It would be of interest to 
investigate the question of the completeness of our solution. In the case of the Heisenberg 
chain, it is known that the analogous ansatz gives complete description of all two-magnon 
states. At this time, it is not clear how to generalize this ansatz for the case of three or 
more electrons, and to find higher conserved currents, but our results give clear indications 
for the integrability of elliptic families of the ID Hubbard models. 



-<Mo [CM + (M+r + r + 2C(A) 

+ip a(r) [t(j - k)e l ^ +p ^ k + t(k - jy^+^'J - CM + r - f - 2((r) 

= E<f> jk - g \t(j - k)e^ + ^ k + t(k - j) e ^+^l , 




(24) 
(25) 
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